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We propose an implementation for quantum logic and com- 
puting using trapped atomic spins of two different species, 
interacting via direct magnetic spin-spin interaction. In this 
scheme, the spins (electronic or nuclear) of distantly spaced 
trapped neutral atoms serve as the qubit arrays for quantum 
information processing and storage, and the controlled inter- 
action between two spins, as required for universal quantum 
computing, is implemented in a three step process that in- 
volves state swapping with a movable auxiliary spin. 
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The field of quantum computing has advanced remark- 
ably in the few years since Shor |IJ presented his quantum 
algorithm for efficient prime factorization of very large 
numbers, potentially providing an exponential speed up 
over the fastest known classical algorithm. Because much 
of today's crytpography Q relies on the presumed diffi- 
culty of factoring large numbers, Shor's discovery has im- 
portant implications to data encryption technology and 
has stimulated much work in the field of quantum infor- 
mation. 

Motivated by this and other theoretical developments, 
there is much interest in identifying and realizing exper- 
imental systems capable of generating large-scale quan- 
tum entanglement. In atomic systems, there have been 
several recent proposals using trapped ions/atoms and 
cavity QED systems f|-[7j]. Indeed, atomic systems capa- 
ble of entangling two qubits have already been realized in 
some of these systems ^,|) . A common element in most 
of these proposals is that the qubits are stored in distin- 
guishably trapped atoms/ions. The proposals differ prin- 
cipally in the nature of the atom-atom interaction (cither 
phonons, photons, collisional, and induced electric-dipole 
moments) and in the way that these interactions are con- 
trolled. 

In this paper, we propose an implementation of a quan- 
tum logic scheme utilizing the direct magnetic spin-spin 
interaction between individually trapped neutral atoms. 
The qubits of this system are stored in the long-lived 
hyperfine ground states of atoms, and coherent control 
of the spin-spin interactions is accomplished by control- 
ling inter-atomic spacings. Our proposal is distinctive in 
that (1) the magnetic spin-spin interaction used to cre- 
ate the inter-atom entanglement is virtually decoherence- 
free and (2) atom-atom interactions are mediated via a 
movable 'header' atom which serves to transfer quantum 
information from one qubit to another (see Fig. 1). The 



header atom can in fact be a different species, and hence, 
in contrast to |5|,|7|,|l0| , the atom trapping potentials are 
not required to be spin-dependent in order to maintain 
trap distinguishability for small atom separations and 
can be realized with far-detuned laser beams. This lat- 
ter distinction is important because near-resonant laser 
traps are a significant source of decoherence. 

We begin by considering in detail the inter-atomic po- 
tential between two neutral atoms separated by an inter- 
nuclear distance R. For the moment, we assume two spin 
1/2 alkali atoms and momentarily neglect the hyperfine 
interactions. The potential can be written as the sum of 
three terms p"T| . 

V(R) = V T {R)V T + V s (R)Vs + V D , 

Vt and Vs are the (electronic) spin triplet and singlet 
potentials, respectively. Vd represents the long range di- 
rect magnetic dipole interaction between two atoms. Vt 
and Vs are the projection operators into the total elec- 
tronic subspace 1 (triplet) and (singlet). The difference 
between Vr and Vs represents the exchange interaction, 
which is typically most important when R is less than the 
LeRoy radius Ro (< 40ao for two identical alkali atoms). 
For two different atom species, the exchange interaction 
is considerably suppressed beyond the contact limit (a 
few ao). In the long range limit both Vt and Vs are 
dominated by the van de Waals term —Ce/R e |0- 

At low energies, we can re-express the first two terms 
of the potential by writing the spin triplet and singlet 
potentials in terms of the scattering lengths, ot and as 



V U (R) = 



a v 5{R), 



v = T,S. 



and explicitly evaluating the projection operators to yield 



V(A) = ^(~*r + ±a 8 )I8(A) 

AirTi 2 . 1 _ _ -* 

H — J7~( a T - a s)^i ■ &2 o{R) + V D . 



(1) 



where a v is electron Pauli spin operators ]T3| . 

At this juncture, we point out that the recent Inns- 
bruck proposal |Hj employs the close-range part of the 
potential represented in the first line of Eq. (|l|) in a type 
of "controlled collision". In our scheme we will use the 
spin-dependent interaction, i.e. the long range atomic 
magnetic interaction represented in the last term of Eq. 
@). The proposal of Brennen, et. al. @ relies on the 



near-resonant electric dipole interaction (not present here 
in the ground state Hamiltonian) . 

It is convenient to re-express the second term of Eq. 
(|I|) by assuming that the two interacting atoms (denoted 
by subscripts q and h) are harmonically bound in cylin- 
drically symmetric traps with characteristic radial and 



axial sizes: a a 



a qz , and a^ r and furthermore that 



the atoms occupy the ground states of their respective 
traps |0) = |0)q|0)/j. In this case, we obtain 



Je = 



Airh z 
' M 



(a T - a s )8(r q - r h - z z)\0) 
a 2 hu! -- 4- 



at a z 



for a reference harmonic trap frequency Huj with ground 
state size a. We have used f q and for the nuclear 
coordinates of the atoms with respect to their own trap 
centers, which are displaced by d = (0, 0, zq), and we have 
defined a v — J a 2 u + a 2 it/ , (y — r,z). It's important to 

recognize that Je decays exponentially with the nominal 
atom-atom separation, zq. 

The last term in Eq. (|l|), Vd, contains three sepa- 
rate terms corresponding to electron-electron, electron- 
nuclear, and nuclear-nuclear magnetic dipole interac- 
tions. Between alkali atoms, the strongest is the electron 
dipole interaction 



T/ee f^e 



[o q ■ °h - 3(i? ■ S q ){a h ■ R)}, 



where fj, e is the Bohr magneton. The strength of this 
interaction is 



le{R) 



R? 



5X10 11 (f)V), 



while 7 en (i?) (electron-nuclear) and j n (R) (nuclear- 
nuclear) are about 10 -3 and 10~ 6 times smaller respec- 
tively. Therefore one may effectively write the spin- 
dependent interaction Hamiltonian as 

H = Je(zq) o q ■ a h + j e (R)[a q ■ a h - 3{R ■ a q ){a h ■ R)]. 

Typically, we will have j e (R) > Je(zq) for R > lOOOao 
between two identical atoms. 

We will now discuss how this interaction Hamiltonian 
can be used for logic gates. We first point out that 
this interaction resembles the quantum gate implemen- 
tation using the Heisenberg spin (exchange) interaction 
||l5f Hj = J{t)a\ ■ (?2, which is known to be universal. 
For Jq dt{J/K) = 7r/4(mod27r), its unitary evolution op- 
erator creates a swap gate 



U, 



S wap(T)|i)i|j)2 = exp(-j-)|j)i|i) 2 . 



which in turn can be use to generate XOR (controlled- 
NOT) gates by incorporating single bit operations 



|15| . However, our interaction Hamiltonian includes an 
anisotropic term. Fortunately, we can borrow a de- 
coupling technique developed in NMR |l6| to effect the 
conversion of a q ■ 3^ — ► On Z ■ ai lz , which is also universal. 
In fact, the phase gate |15| in terms of these operators is 
simply 

[/ phasc = e i ( 7r / 4 ) cr i^2; x giW 4 )^ x e i(7r/4)<T 2i ^ 

from which Uxor can be easily made [plJl^]. Further- 
more, the swap gate, which we will require, can be made 
according to 

C/swa P (1^2) = C/xor(1,2)C/xor(2, 1)C/ XO r(1, 2), 

where Uxon{i,j) denotes a C-NOT with i as the control 
bit operating on j. The necessary de-coupling is achieved 
through a "stirring" radio frequency field acting only on 
the /i-atom [|l6| , and is most easily discussed in the con- 
text of the following model Hamiltonian 

H s (t)= riuj x a lz + hu 2 a 2z + n s (a 2+ e- iuJst + h.c.) 

+7e(i?)^l ■ ( ?2-3 ( Ji z( 7 2z (i?-z) 2 ], (2) 

where ujs is the frequency of the stirring field, and £ls is 
the Rabi frequency of the stirring field. By analyzing this 
system in the rotating frame defined by Un — e lwst<T2z , 
we obtain jl6| U^u 2 ±U-n — > a 2+ e ±l( ^ Lt , and invoking a 
rotating wave approximation, the desired result is ob- 
tained, 



s ~ 7eOR)[l - 3(i? • z) 2 ]a lz a 2z 

+ kuj 1 a lz + h(u 2 - UJs)&2z + ^s(ct 2 4 



cr 2 _. 



(3) 

Although there are unwanted single atom terms in the 
second line of this Hamiltonian, they can be easily com- 
pensated with one-bit rotations. For our system, a simi- 
lar procedure yields the following effective Hamiltonian, 



n cS {R) 



[Je(z q ) + j e (R) 

Je(zq) O qz Ohz- 



3 7e (i?)(i? ■ z) 2 ]a qz a h 



(4) 



Interestingly, we note that the spin and spatial depen- 
dence of the operators factorizes. This implies that co- 
herent spin-spin interactions only require that the mo- 
tional states of the atoms remain unchanged — the atoms 
are not necessarily required to be in the ground state 
|0) of their respective trapping potential. This particular 
feature of our proposal will be discussed in detail else- 
where. At zq > Rq, the effective spin-spin interaction 
strength is 

Mzo)~(je(R)[l-Hz-R) 2 }), 

which for atoms in the ground states |0) previously de- 
scribed is readily evaluated 



1 



3(i-i?) 2 ] 
1 



v / 2^a 2 z 2af 

(a 2 + z 2 ) 



dz exp 



2a 2 



2z 



■ exp 



z 

2aJ 



erfc 



\72a r 



where erfc(.) is the complementary error function. The 
geometry of the system of two interacting spins are il- 
lustrated in Fig. g. The result of the effective interac- 
tion is shown in Figure]^, we note that Je(zo) is in the 
kHz range for a distance of lOOOao (~ 50 nm), which 
will be more than adequate for gate operations for atoms 
trapped in far off-resonant optical lattices. 

The principle challenge in implementing this scheme 
is in providing the appropriate confining potentials for 
the atoms. On one hand, the trapping potentials for the 
individual atoms need always be distinguishable in order 
to maintain identifiable qubits. On the other hand, as 
we can see from Figure 2, the atoms need to be in close 
proximity (~ 50 nm) in order for an appreciable inter- 
action rate even for this 'long-range' potential. Previous 
proposals also requiring small inter-atomic spacings have 
suggested spin-dependent traps created by optical lat- 
tices with polarization gradients (|,|7],|l0| . Because of the 
nature of these types of traps, the types of atomic ma- 
nipulations are rather restricted, and hence scalability is 
difficult. 

To circumvent this complication, we will use two dif- 
ferent atomic species, one for the (stationary) quantum 
register, and one for the quantum header atom. Each 
species of atom will be separately trapped by different 
laser fields. By appropriate choice of atom and frequency 
of the trapping fields, we can make these traps essentially 
independent. For a concrete example, consider a quan- 
tum register consisting of an array of single atoms (type q 
for qubit) trapped in 3D standing wave formed by inter- 
ferin g la ser field of a CO2 lasers (wavelength Aco 2 ~ 10-6 
fim) ]18| . The qubits will be separated by Aco 2 /2 which 
is more than enough to allow individual addressing, and 
at this separation, the long range Casimir-Polder inter- 
action is negligible |l2|]. The trapping details are dis- 
cussed in the appendix, but we point out that the poten- 
tial V is very well approximated by the dc-polarizability 
of the atom a (0)and the laser electric field amplitude E 
as V = —a(0)E 2 /2. The relevant parameters are tabu- 
lated for alkali atoms in Table |. 

A separate laser field provides confinement for the 
header atom (of a different type atom, h). By choosing a 
trapping wavelength somewhat closer to the atomic res- 
onance of h (and detuned to the blue of the resonance), 
we can provide a potential which acts principally on the 
h atom. The potential depth for this trap is given by 
Knax = ^l/4<5l, with f2 the Rabi frequency of the 
laser, and Sl — ^>l~ ^0 the detuning. The h atom is also 
affected by the far off-resonant CO2 laser field of course, 
but we can arrange for the off-resonant blue detuned field 
to dominate the far off-resonant CO2 laser potential by 
suitable choice of atoms. Trapping parameters for this 
case are listed in Table Q. The quantum register atoms 
(type q) will also be affected, at some level, by the blue 
lattice, but the detuning between the blue field and the q 
atoms will be much larger, so the potential will be dom- 
inated by the CO2 laser field for the q atoms. 

Gate operations in this system can be achieved in a 



three step process requiring quantum state swapping be- 
tween the quantum bits and the header atom. To exe- 
cute a gate operation between two qubits and qj, we 
first translate the header atom h to the location of 
and perform a state swap q^h. The header atom is then 
translated to site qj and the gate operation between h(qi) 
and qj is performed. Finally the header atom is trans- 
lated back to qt to and the state swap is repeated. The 
header atom effectively acts as an quantum bus between 
the qubits, and in this sense our scheme shares certain 
features with the quantum gear machine proposed by Di- 
Vincenzo E3 . 

Single-bit operations can be realized either by directly 
addressing the individual qubits or, alternatively, we 
can use the header atom as a mediator. The latter option 
may be easier in some cases than the direct spatial selec- 
tion of qi because the h atoms can be sparsely distributed 
and have different resonance level structures. The sin- 
gle bit operation will again be a three step process: 1) 
perform a state swap between qi and h, 2) perform the 
arbitrary qubit operation on h, 3) repeat the state swap 
between h and q t . 

In considering the ultimate scalability of this, and 
other lattice-based schemes, it is necessary to compare 
the characteristic intrinsic decoherence time of the sys- 
tem to the gate time plus the transport time of the mov- 
ing atoms [£0| . Additionally, the transport of the moving 
atoms (the h-type atom in this case) must be adiabatic 
such the motional state of the atom remains unchanged. 
This latter condition implies constraints on the magni- 
tude of the motion, which we can estimate using pertur- 
bation theory. Consider the Hamiltonian for a one di- 
menstional harmonically trapped particle with mass M 
and trap frequency u>t subjected to a force F(t), 

H=^ + hdu 2 q 2 -qF(t), (5) 

adiabaticity condition for the header qubit translation re- 
quires its motional state wave-function to be essentially 
unchanged. This problem is equivalent to the problem 
of a translating SHO potential Moj 2 [q — q (t)] 2 /2 [with 
F(t) — Mu>fqo(t)] up to a deterministic phase factor due 
to Mu) 2 q 2 t (t)j2. Calculating the probability for excita- 
tion out of the ground state is a standard textbook prob- 
lem and the result to first order is 

W = 2 J x exp(- Wt r), (6) 

for a time dependent force F(t) — (F t/lj)/(t 2 + t 2 ). 
^M(Sv) 2 is the energy gained from the impulse M5v — 
F(t)dt of the force. We see that adiabaticity is main- 
tained even after the translating atom gains a very large 
speed, but satisfying the condition u>tT » 1, i.e. a force 
to be slowly turning on and off compared with the har- 
monic trap period. Similar conclusions are reached for 
an initial coherent motional state wave-packet. This con- 
dition effectively then puts no constraint on the header 



atom speed, contrary to the strong conditions as obtained 
in Ref. ^0|. For our problem, creative pulse shape de- 
sign will allow the header atom to be adiabatically trans- 
ported over many qubits within the single photon scat- 
tering coherence time. 

Our discussion thus far has been limited to alkali atoms 
with no nuclear spin (e.g. 78 Rb). When the nuclear spin 
I is nonzero, the atomic spin takes on values F = I ± 
1/2 and we must include the hyperfine interaction Vhf ~ 
a\ l i<j-d n . The spin- spin interaction becomes considerably 
richer in detail. However, if a strong Zccman interaction 
is applied using a uniform magnetic field, the resulting 
two manifolds of Zeeman states correspond roughly to 
the electronic spin up/down such that the good basis 
becomes \I,S,I Z ,S Z ) ©J. Alternatively, we could chose 
an atom with no nuclear spin such as Rb (radioactive 
lifetime about 20 minutes). 

In summary we have proposed a new quantum comput- 
ing implementation with trapped atomic spins. Utilizing 
dual optical lattices for two different type of atoms pro- 
vides a novel method to control the binary interaction 
between any pair of qubits. In addition, our proposal, 
being based on the periodic structure of optical lattices, 
is readily scalable, and in particular, redundant paral- 
lel processing of information can be implemented using 
multiple h-type atoms operating on repetitive blocks of 
g-type atoms. This may be useful in implementing error 
correction pal , concatenated coding, and fault tolerant 
computing p3| . 
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topes. This work is supported by the ARO/NSA grant 
DAA55-98-1-0370 and by the ONR research grant No. 
14-97-1-0633. 



«osc («o) 


778 


573 


433 


347 


295 


Ao (nm) 


670 


589 


766 


780 


852 


E R (kHz) 


64 


25 


8.7 


3.7 


2 


no 


0.39 


0.32 


0.19 


0.15 


0.11 


Vco 2 


0.025 


0.018 


0.014 


0.011 


0.009 


TABLE 


II. 


Parameters 


for 


a 


'blue' lattice with Ql 


- 1.6 x 


10 10 (Hz) 


(~ laser 


power of 


10 kw/cm 2 ) 


, 7 = 10 7 (Hz), and S L = 2 x 


10 12 (Hz) 


. For the 



near resonant 'blue' lattice on the ft— type atoms, the effective 

2 f2 2 

single photon scattering rate is approximately j e f[ = r\ 753-7. 



We see as indicated in Table |H] the confining frequency iw is 
indeed much larger than that of CO2 laser (on g-type atoms). 





Li 


Na 


K 


Rb 


Cs 


M 


6.9 


23 


39 


87 


133 


v osc (kHz) 


4061 


2530 


1494 


982 


727 


Hose (ao) 


254 


176 


176 


145 


137 


V 


0.13 


0.1 


0.076 


0.06 


0.05 


Tefi (Hz) 


2.5 


1.6 


0.9 


0.6 


0.5 



APPENDIX A: POSSIBLE TRAP PARAMETERS 



TABLE I. Parameters for different alkali-metal atoms in- 
side a CO2 lattice with intensity / = 10 6 (watts/cm 2 ). 
For the 'red' CO2 lattice, the maximum level shift is 
V max = a(0)^o/4. At an intensity of ~ 10 6 (watts/cm 2 ), 
the single photon scattering rate can provide decoherence 
times of many minutes. Assuming a harmonic approxima- 
tion, the oscillation frequency u OBC inside the CO2 trap is 
^osc = 2^0^ with E%° 2 = h/(2M\ 2 C02 ) the recoil 
energy (in Hz) for emitting/absorbing a C02-photon. The 



Lamb-Dicke parameter 770 = koa osc — y/ E R /v osc (?7co 2 ) mea- 
sures the trap ground state size a osc in terms of the resonant 
wavelength Ao (CO2 laser Aco 2 )- 





Li 


Na 


K 


Rb 


Cs 


M 


6.9 


23 


39 


87 


133 


a(0) (og) 


159.2 


162 


292.8 


319.2 


402.2 


V max (MHz) 


181 


185 


334 


364 


458 


v osc (kHz) 


432 


239 


247 


172 


156 



Computation, (World Scientific, London, 1998), edited by 
H.-K. Lo, S. Popes cu, and T. P. Spiller; A. Yu. Kitaev, 
( |quant-ph/970702l| ). 
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FIG. 1. A one dimensional illustration. The q— type atom 
array are trapped in a red periodic CO2 laser lattice. The 
movable header atom is trapped in a blue lattice. 




FIG. 2. The geometry of interacting header atom and qubit 
atom pair. The large ellipses denote trap ground states with 
trap centers crossed and separated by zq. Solid circles with 
arrow heads denote electron spins separated by R. 
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FIG. 3. The solid line denotes the exchange interaction Je 
assuming a absolute difference of the \as — ar\ = 100 (ao), 
while the dots are are numerical results of the averaged spin 
dipole interaction strength Je for a qr — a qz — 400ao and 
cihr = a,hz = 100ao. The dotted line represents the simple 
1/zq dependence. As pointed out in the text, the exchange 
interaction will be significantly suppressed for the case where 
the two atoms are two different species. 



